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Marginal  Stability  of  Ion-Acoustic  Waves  in  a 
Weakly  Collisional  Two-Temperature 
Plasma  Without  a  Current 


1.  INTRODUCTION 


For  collisionless  current-free  and  current-carrying  two-temperature  plas¬ 
mas,  the  stability  of  ion-acoustic  waves  has  been  studied,  and  the  results  have 
been  accepted  for  many  years.  For  example,  linear  Vlasov  theory  predicts  that 


ion-acoustic  waves  in  current -free  plasma  are  Landau-damped,  and  the  plasma  is 


1 


stable.  In  the  presence  of  a  current,  ion-acoustic  waves  become  unstable  when 
the  current  exceeds  a  threshold  value.  * 


For  collisional  plasmas,  the  stability  situation  is  unclear.  On  the  one  hand, 

there  are  studies  of  collisional,  current-free  plasma  that  indicate  that  ion-acoustic 

2  3  4  2  3 

waves  are  damped.  ’  '  In  their  work,  Kulsrud  and  Shen  and  Ono  and  Kulsrud 


assumed  that  the  electrons  were  isothermal  and  solved  the  Fokker-Planck  equation 


.4 


for  the  ions  to  find  that  i-i  collisions  damped  ion-acoustic  waves.  Buti  solved  the 


(Received  for  publication  5  August  1987) 

1.  Fried,  B.D.  ,  and  Gould,  R.W.  (1961)  Longitudinal  ion  oscillations  in  a  hot 
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plasma,  Phys.  Fluids,  _4_:139-147. 

Kulsrud,  R.M.,  and  Shen,  C.S.  (1966)  Effect  of  weak  collisions  on  ion  waves, 
Phys.  Fluids,  9:177-186. 

Ono,  VI. ,  and  Kulsrud,  R.M.  ( 1975)  Frequency  and  damping  of  ion  acoustic 
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Fokker-Planck  equation  in  the  weakly  collisional  limit  for  both  the  electrons  and 
the  ions  to  find  that  e-i  collisional  undamping  was  small  compared  with  i-i  colli¬ 
sional  damping  for  a  wide  range  of  wavenumbers  and  temperatures.  On  the  other 
hand,  there  are  studies  of  weakly  collisional  plasma,  driven  unstable  by  the  pres¬ 
ence  of  a  current,  which  indicate  that  e-i  collisions  have  a  significant  undamping 
2  5 

effect  '  on  ion-acoustic  waves  and  lower  the  critical  current  required  to  produce 


a  current-driven  instability.  In  these  studies,  however,  the  authors  gave  no 
discussion  of  the  stability  of  the  plasma  in  the  limit  of  zero  current.  In  our  view, 
the  stability  of  ion-acoustic  waves  in  a  weakly  collisional,  current-free  plasma  is 
unresolved.  The  purpose  of  this  report  is  to  present  a  resolution  based  upon  the 
Balescu-Lenard -Poisson  equations. 

The  Balescu-Lenard  kinetic  equations  for  the  one-particle  distribution  func¬ 
tions  for  a  weakly  coupled  (g  <<1)  plasma,  which  may  be  used  to  study  long  wave¬ 
length  phenomena,  are  derived  by  truncating  the  Bogoliubov-Born-Green-Kirkwood- 
Yvon  hierarchy  to  the  first  order  in  g,  by  applying  the  adiabatic  approximation  and 

g 

by  using  time-asymptotic  solutions  for  the  pair  correlation  functions.  Here,  g  is 
the  plasma  parameter.  Use  of  the  adiabatic  approximation  implies  that  the  pair 
correlation  functions  relax  on  a  time  scale  that  is  fast  compared  with  that  of  the 
one-particle  distribution  functions,  and  use  of  time-asymptotic  solutions  for  the 
pair  correlation  functions  implies  that  plasma  waves  are  not  growing  in  time.  Al¬ 
though  the  adiabatic  approximation  is  known  to  break  down  for  to  greater  than  the 
electron  plasma  frequency  (Jasperse  and  Basu,  Section  IVC),  we  assume  that  it 
is  valid  for  low  frequencies  and  long  wavelengths.  We  also  assume  that  the  time- 
asymptotic  solutions  for  the  pair  correlation  functions  are  valid  in  stable  plasmas 
up  to  the  limit  of  marginal  stability. 

IrL'this  report,  wef' present. two  important  results  for  a  weakly  collisional, 
current-free,  two-temperature  electron-ion  plasma:  (1)  that  e-i  collisions  have 
an  undamping  effect  on  ion-acoustic  waves,  and  (2)  that,  for  appropriate  ^values  of 
the  plasma  parameters,  ion-acoustic  waves  become  marginally  stable.  (Our  re¬ 
sults  are  based  on  the  closed-form  solution  for  the  dielectric  function  for  the  line¬ 
arized  Balescu-Lenard-Poisson  kinetic  equations,  which  is  given  by  Eq.  (2).  The 
weakly  collisional  ordering'we  consider,  which^is  necessary  to  obtain  an  iterative 


5.  Ong,  R.S.B.,  and  Yu,  M.Y.  (1969)  The  effect  of  weak  collisions  on  icn- 

acoustic  wave  instabilities  in  a  current-carrying  plasma,  J.  Plasma  Phys. 
JL425-433;  Stefant,  R.  J.  (1971)  Influence  of  electron-ion  collisions  on  ion 
acoustic  waves,  Phys.  Fluids,  14:2245-2246. 

6.  Montgomery,  D.  C. ,  and  Tidman,  D.A.  (1964)  Plasma  Kinetic  Theory, 

McGraw-Hill,  New  York.  "  ~ 

7.  Jasperse,  J.  R.,  and  Basu,  B.  ( 19G6)  The  dielectric  function  for  the  Balescu- 

Lenard-Poisson  kinetic  equations,  Phys.  Fluids,  29:110-121. 
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solution  lor  the  collisional  dielectric  function,  imposes  two  independent  restric¬ 
tions  on  the  parameters  that  characterize  the  plasma.  "They  are:  (1)  u  g./w  .  = 

(  v  .  /  to  .)  (k  /k)  <  1  or  1  /kX  .  =  v  .  /kv™.  <  (m  /m.)*  2;  and  (2)  v  .. /  to .  <<  1  or 

ei  pi  e  ei  ei  i  e  e  i  n  la 

1/kA  ..  =  v../kv,T,.  <  1.  Here,  v  .(v..)  is  the  e-i  (i-i)  collision  frequency;  w.  is 
n  ii  Ti  ei  u  n  J  la 

the  ion-acoustic  wave  frequency;  to  ^  is  the  ion  plasma  frequency;  kg  is  the  electron 
Debye  wavenumber;  Xg.(X..)  is  the  e-i  (i-i)  mean  free  path;  vq’e(V'pj)  is  the  elec¬ 
tron  (ion)  thermal  speed;  and  the  other  quantities  have  their  usual  meanings.  The 
approximate  marginal  stability  condition  is 

V“ia  ’<  W  Vkl  '  (2")1/2 


+  <Te0/T10»3/2  exp(-3/2-Te0/2Ti0H. 


(1) 


provided  that  T-^/T^  and  k/kg  <<1  and  the  inequality  1  >  vei/a)ia  holds.  Here, 
1/2  3/2 

vei  =  (2  In  A0i/ 12 tt  )  g  w  #  where  g  is  the  plasma  parameter.  In  is  the 

Coulomb  logarithm,  and  w  is  the  electron  plasma  frequency. 


2.  CLOSED-FORM  SOLUTION  FOR  THE  COLLISIONAL 
DIELECTRIC  FUNCTION 

We  begin  with  the  nonlinear  Balescu-Lenard -Poisson  kinetic  equations  for  an 
unmagnetized,  two -temperature  electron-ion  plasma.  We  consider  a  quasi-steady, 
zero-order  state  where  the  zero -order  distribution  functions  are  Maxwellian  with 
unequal  temperatures  (Tg0  /  T.Q).  In  order  to  obtain  the  dielectric  function  that 
contains  the  effects  of  e-e,  e-i,  i-e,  and  i-i  collisions  to  the  first  order,  we  line¬ 
arize  the  equations  for  small  amplitude,  electrostatic  perturbations  about  the 

O 

zero-order  state  and  apply  the  collisional  propagator  expansion  method  appropri¬ 
ately  generalized  for  a  two -constituent  plasma.  Using  the  techniques  discussed  in 
7 

Jasperse  and  Basu,  we  find  that  all  of  the  integrals  that  appear  in  the  dielectric 
function  may  be  reduced  to  the  simpler  forms  presented  here.  Since  the  algebra 
involved  is  straightforward  but  lengthy,  the  details  of  the  calculation  will  be  re¬ 
ported  elsewhere.  The  dielectric  function  to  the  first  order  in  collisionality  is: 

c(1)  (k.  W)  =  1  +1  (k  /k)2  {W0(za)  +  i  Ina6  [AaB(za)  +  Ba6(za)]}.  (2) 

a  g 

8.  Jasperse,  J.  R.  (1984)  A  propagator  expansion  method  for  solving  linearized 
plasma  kinetic  equations  with  collisions,  Phys.  Letters,  106A:379-382. 
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9  2  2 

Here,  k“  =  4tt  nAq  /T  ,  k  =  Ik  I,  z  =  w/kv™  ,  v™  =  T  /m  ,  nryp 
a  CP  a  aO  1  a  Ta  Ta  aO  a  nb 


Va6  /kvT a*  vaB  =  (4/3)(2l)1/2  < Wo^'f  TS>  (1 


4  ™  M>(1  4  C«P)_ 


2  2 

i  n  is  the  a  -  6  collision  frequency,  caR  =  v^,p  /vTa  ,  =  m(i  /m,.  ,  and 

AaPU>  =  (|HH  =aB  ^(l+m^r1  ||l+3ea6  (T^  -  1HH  'af>  I-', 

-  !<“  +  *  ,<»(1  +E«B  1-1 

-  <r,)  +  l(T„e  -  Dj”!  <011. 

B«P  f'->’3c«e  (TaS  -»«  "I  *,*.,,<>  <>♦'«»  )3/2(1  +  m,.|i,'ll“)(C». 

1 

1{"1  (c)  =  -c1_n  /0dM  u2(l  -m2+  ectB)"(4'n)/2W2.n  la^(u)c], 

J(n,)  (  o  =  -  r,  J  du  ti2(l  -  u2  +  caP  )“n  f  dt  t1  "  n 
aB  J0  0 

x  exp  [ ib  ^  (  u)  c  t  -  c^p ,  (  u  )t2]  W'n  [a  ^  (u  )  C  +  i e  ( ...  w2t], 

2  I/O 

In  the  above  n  =  0,  1,  T  .  T  ag/TpQ,  a^p  =  1(1  ■<  £,,)/(!  “  11  “  p  ,,.•)] 

q  .  =  q  ./qrt»  b  .  =  |(1  •>  ,..)d  -  x  ^  )]]  /2,  and  cap  =  (^)  eap  x 

[(1  -  u“)(  1  -  c^p  )"  +  £  1 1,  id].  The  W^-functions  are  given  by 

OO 

Wn<  C)  =  in  jQ  dt  tn  '  1  exp  (i-:t  -  t2/2), 

where  n  =  0,  1  ,  2,  .  .  . ,  VV  ( r )  -  (<l/d  )n  W  (  r,),  and  \V  (  r)  is  related  to  the  plasma 

n  u  1/9  "i  /  9 

dispersion  function  /  (  ’. )  by  \\  (  r )  1  “  ('  /2  ~)  Z('/2  “).  In  Kq.  (2),  k_  has  been 
oriented  along  the  positive  z-axis.  We  note  that,  when  T0q  r  T.^  in  Llq.  (2),  we 
obtain  Kqs.  (1)  and  (2)  of  Jasperse  and  Hasu.  ^  Also,  when  T  ^  =  T.^  and  ion  mo¬ 
tion  is  neglected,  we  obtain  Eqs.  (26)  through  (20)  of  Jasperse  and  Hasu*  for  the 
high-frequencv  dielectric  function. 


f.  Jasperse,  J.H.,  and  Hasu,  11.  (1087)  (  ollisional  enhancement  of  low -frequence 
densitv  fluctuations  in  a  ueaklv  (  ollisional  electron-inn  plasma,  l’hys.  Rev. 

1 . ette r s ,  58:1422-1425. 
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3.  STABILITY  ANALYSIS 

Numerical  solutions  of  the  dispersion  relation,  £  ^(k,  w)  =  0,  for  ion-acoustic 
waves,  that  is,  for  Vrp.  <  | oj/k j  <  v,pe  and  k/kg  <1,  and  for  Aag  =  12  n  /g,  are 
presented  in  Figures  1  and  2,  Figure  1  shows  the  dependence  of  Lm  o>  /Re  u>  (Im  w 
<  0  means  stability)  on  k/kg  and  T^/TVq  for  an  Argon  plasma  with  g  =  2  x  10“®. 


.20  .2! 


Figure  1.  Values  of  Imw/Re  w  From  a  Numerical  Solution  of 
the  Dispersion  Relation  for  Ion-Acoustic  Waves  in  a  Weakly 
Collisional  Argon  Plasma  Using  Eq,  (2) 
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In  what  follows,  we  present  some  analytical  results  that  we  obtain  from 
Eq.  (2)  by  making  approt  riate  approximations.  From  Eq.  (2),  we  see  that 
contains  eight  collisional  terms,  two  terms  for  each  of  the  four  collisional  pairs: 
e-e,  e-i,  i-e,  i-i.  The  imaginary  parts  of  the  terms  for  each  of  the  four  colli¬ 
sional  pairs,  denoted  by  x  ^ g  ,  contribute  to  the  imaginary  part  of  the  wave  fre¬ 
quency  as  damping  or  undamping.  For  ion-acoustic  waves  at  long  wavelengths,  we 
find  that  X  ^  :  X  ^  :  X  ^  :  X  ~  1  *  (16/5)  (m^m.)1  /2  (T^/T iQ) 1  /2 :  me/mi; 


1 :  (16/5)  (m  /2m . ) 

e  i 


m  /m..  We  also  find  that  (m  /m.)  J  . 
e  1  e  i  ei 


,  (m  /m.)  I  J(1)  |  and  e  .  j  I(1>  |  «.'I(0). 
*  e  i  ei  ei  '  ei  1  ei 


Hence,  for  m  /m.  <<  1  and  (T  n/T  ..)^2  <  <  (5/16)  (2m .  /m  )*^2,  £  ^  may  be  well 

O  X  6U  lU  1  0 

approximated  by 


e(1)(k.  u)  )  =  1  +  l  (k  J k)2  W  (z0)  4  (i/2)  (ke/k)2  pei  4°/  (z  >. 
a 


In  addition,  we  find  that  for  <  <  |u)/k  |  <<  v™  ,  I  .  ~  -  1/z  ,  and  that  e  *  may 

x  i  i  e  ei  e 

be  approximated  asymptotically  by  using  the  appropriate  expansions  for  W  : 

e(1)  (k,  to )  ~  1  +  (k  /k)2  -  (to  ./to)2  (1  +  3z72)  +  i(k  /k)2 
*’**  6  pi  i  e 


x  { (  u/2)1  (m./m  )A/  z  [(m  /m.) 

l  e  e  e  1 


+  (TeO/TiO)3/2  exP( (^ )]  -  (ve./2w) 


Ion-ion  collisions  result  in  a  damping  term  in  the  dielectric  function  that  may  be 

2  -4 

approximated  analytically  by  (8i/5)  (T^/T-g)  (v../ to)  (kg/k)  z  .  and  can  be  ne¬ 
glected  compared  to  the  e-i  undamping  term  for  (T  /T  )^^  <<  (5/16) 

1/2  e  u  i  u 

(2rm/me)  .  Solving  the  dispersion  relation  for  k/kg  <<  1  and  |lmto|<<  Re  to  us¬ 
ing  Eq.  (3),  we  obtain 


Re-  >kHT(0OT10l/m.|''<=V 


Im  to  -  -  a) .  {(tt/8)^2  [(m  /m.)^2 

la  e  l 


+  <TeO/TiO>  exP("3/2  -TeO/2TiO)]  “  (  vei/4“ia)  i  » 


which  are  valid  when  T...  <<  T  n  and  v  .  <w .  .  In  Eq.  (4),  we  have  neglected  the 
first-order  collisional  shift  in  the  expression  for  Re  to.  Eq.  (5)  for  Im  to  is  the  sum 
of  three  terms:  The  first  is  electron  Landau  damping,  the  second  is  ion  Landau 


flares 


damping,  and  the  third  is  an  undamping  term  due  to  e-i  collisions.  From  Eq.  (5), 
we  find  the  marginal  stability  condition  given  in  Eq.  (1).  When  T  q/T.q  >  20  and 
k/ke  <  0.  1,  Eq.  (5)  yields  values  that  are  within  10  percent  of  the  results  shown  in 
Figures  1  and  2.  For  smaller  values  of  Teo/TiO'  numerical  solutions  are  neces¬ 
sary  for  accurate  results. 

By  taking  the  velocity  moments  of  the  perturbed  electron  distribution  function 
and  neglecting  collisionless  damping  effects,  we  find  that 

Sel/n0  =  ‘  (i/2)  (vei/u))J  (e^l/TeO)’ 

Ll/Lo",i/2)Wel/",<eVTeO>- 

where  n  j,  tjL,  and  Tgl  are,  respectively,  the  perturbed  electron  density,  per¬ 
turbed  potential  and  perturbed  electron  temperature  in_k_-  to  space.  These  results 
show  that  the  electrons  are  not  isothermal  and  that  there  is  a  collision-induced, 
temperature  perturbation  associated  with  the  ion-acoustic  mode.  We  recall  that 
in  the  linear  Vlasov  theory  of  the  ion-acoustic  mode,  the  electrons  are  isothermal 

and  T  =  0. 
el 

We  have  performed  a  quasilinear  analysis  of  the  Balescu-Lenard -Poisson 
equations.  We  find  that 


(6) 

(7) 


Te0(t)  St 


-  Teo'» 


-  J  dk  dm  u>k)  I  e$j(t)  /Te0(t) 


(8) 


where  Im  to  is  given  by  Eq.  (5).  For  stable  ion-acoustic  waves,  Eq.  (8)  indicates 

K 

that  T  initially  increases  with  time  while  the  electric  field  energy  density, 

^eu  ^  i  2  t  ^  1 2 

k^  cp  j(t)  |  /8  tt,  initially  decreases  since  3  |  tp ^(t)  |  / 3t  =  2  Imw  ^  |  cpj(t)  |  <0.  A 
more  detailed  discussion  will  be  presented  elsewhere. 


4.  CONCI.l  SION 


In  conclusion,  we  have  shown  by  using  the  Balescu-Lenard -Poisson  equations 

that  e-i  collisions  have  an  undamping  effect  on  ion-acoustic  waves  for  weakly  col- 

lisional  plasma.  The  fact  that  this  undamping  effect  persists  even  when  TgQ  =  T 

is  discussed  in  Jasperse  and  Basu.  ®  For  thermal  equilibrium  (TgQ  =  T.q)  plasma, 

the  undamping  effect  of  e-i  collisions  on  thermal  fluctuations  at  the  ion-acoustic 

g 

frequency  from  Jasperse  and  Basu  is  1.2  (ve^/4).  For  nonequilibrium  (Teg>>T  q) 
plasma,  the  undamping  effect  of  e-i  collisions  from  Eq.  (5)  is  '->e./4.  see  that 
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the  undamping  effect  operates  for  equilibrium  as  well  as  nonequilibrium  plasma, 

_3  /9  '  1 

and  has  a  value  which,  except  for  the  TpQ  dependence  of  v  is  nearly  independ¬ 
ent  of  TgQ  and  T.^.  The  difference  between  equilibrium  and  nonequilibrium  plasma 
is  that,  for  equilibrium  plasma,  ion  Landau  damping  easily  overcomes  e-i  colli- 
sional  undamping,  and  the  plasma  is  stable,  whereas  for  nonequilibrium  plasma, 
a  wavelength  can  be  found  where  e-i  collisional  undamping  balances  ion  and  elec¬ 
tron  Landau  damping,  and  the  plasma  is  marginally  stable.  Our  analysis  suggests 
that  when  vei/wia  exceeds  the  quantity  on  the  right-hand  side  of  Eq.  (1)  a  collision- 
driven  ion-acoustic  instability  occurs.  Although  the  Baleseu-Lenard  collision 
operator  may  not  be  valid  for  unstable  waves,  it  is  our  belief  that  such  an  instabil¬ 
ity  exists  but  the  actual  threshold  condition  may  be  slightly  different.  In  order  to 
test  this  idea,  we  suggest  that  an  experiment  be  performed.  The  threshold  param¬ 
eters  given  in  Figure  2  should  be  useful  in  designing  such  an  experiment. 
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